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Abstract 

This document emerged from work that started in January 2012 as a joint effort by the ATLAS, CMS 
and LPCC supersymmetry (SUSY) working groups to compile state-of-the-art cross section predic- 
tions for SUSY particle production at the LHC. We present cross sections for various SUSY processes 
in pp collisions at yfs = 7 TeV, including an estimate of the theoretical uncertainty due to scale varia- 
tion and the parton distribution functions. Further results for higher LHC centre-of-mass energies will 
be collected at https://twiki.cern.ch/twiki/bin/view/LHCPhysics/SUSYCrossSections. For squark and 
gluino production, which dominate the inclusive SUSY cross section, we employ calculations which 
include the resummation of soft gluon emission at next-to-leading logarithmic accuracy, matched to 
next-to-leading order (NLO) SUS Y-QCD. In all other cases we rely on NLO SUS Y-QCD predictions. 



1 Introduction 



The search for supersymmetry (SUSY) is a central activity of the LHC physics programme. In the framework 
of the Minimal Supersymmetric extension of the Standard Model (MSSM) with R-parity conservation, SUSY 
particles are pair produced. At the LHC, the most copiously produced SUSY particles are expected to be the 
strongly interacting partners of quarks, the squarks (q), and the partners of gluons, the gluinos (g). The squark and 
gluino pair-production processes to be considered are 

pp -> qq, qq* , qg, gg + x , (l) 

together with the charge conjugated processes. In Eq. (G3 the chiralities of the squarks, q = (qL, qn), are sup- 
pressed, and we focus on the production of the partners of the (it, d, c, s, b) quarks which we assume to be mass- 
degenerate. The production of the SUSY partners of top quarks, the stops (t), and, when appropriate, the partners 
of bottom quarks, the sbottoms (b), has to be considered separately due to parton distribution function (PDF) ef- 
fects and potentially large mixing affecting the mass splittings. In this case, we explicitly specify the different mass 
states in the pair-production processes, 

pp^tJtMK + X * = 1,2, (2) 

where i = 1, 2 corresponds to the lighter and heavier states, respectively. Furthermore, if weak-scale SUSY is 
realised in nature, the production of SUSY partners of the electroweak particles in the Standard Model, namely the 
gauginos and higgsinos, need also to be considered. 

A comprehensive search program sensitive to the production of SUSY particles has been underway at the LHC, 
since the beginning of the data-taking process. These searches aim for a broad range of possible final states [Q]|2]. 
Given the importance of SUSY searches at the LHC, accurate knowledge of theoretical predictions for the cross 
sections is required. The interpretation, in terms of exclusion limits, of the results from the Tevatron, the 2010 
LHC data and the first set of the 2011 LHC events, used next-to-leading order (NLO) predictions. However, 
over the past few years, substantial progress has been made in calculating higher-order corrections for squark 
and gluino production cross sections. In particular, the resummed results at the next-to-leading logarithmic (NLL) 
accuracy matched to NLO predictions have become available, allowing for a significant reduction of the theoretical 
uncertainties related to unknown higher perturbative orders. Therefore, it is desirable to implement this recent 
progress in the forthcoming set of analyses of the most recent LHC data sets. In the case of electroweak or strong- 
weak production, i.e. production processes in which sleptons and gauginos are involved, the NLO predictions are 
used. 

This document emerged from discussions among the SUSY working groups of ATLAS, CMS and of the LHC 
physics center at CERN (LPCC) E). The aim is to provide a reference for the collaborations on the evaluation 
of the production cross sections and their associated uncertainties, along the lines of similar recommendations 
developed by the PDF4LHC |4| forum and by the Higgs cross section working group |5]|6]. For the purpose of 
illustration, this document will present some explicit results for y/s — 7 TeV. The detailed cross section values 
for the relevant processes and SUSY models considered by the experiments, as well as the results for higher LHC 
centre-of-mass energies, will be collected at the SUSY cross section working group twiki page 0. 

The next section briefly describes the current state-of-the-art higher order calculations, followed by the prescription 
used for the treatment of theoretical uncertainties in section [3] In sections @] and |5j the production cross sections 
for coloured SUSY states and electroweak SUSY productions are shown, respectively. A summary of the results 
and of the future prospects is given in section [6] 

2 Higher order calculations - NLO+NLL 

The dependence of hadron collider observables on the renormalisation and factorisation scales is an artifact of 
perturbation theory. This is typically reduced as higher-order perturbative contributions are included. Assuming 
that there is no systematic shift of an observable from order to order in perturbation theory, for example due to the 
appearance of new production channels, the range of rates covered by the scale dependence at a given loop order 
should include the true prediction of this rate. The scale dependence gives us therefore a lower limit on the theory 
uncertainty of a QCD prediction, which becomes smaller as higher-order SUSY-QCD corrections are included. To 
estimate the scale uncertainty in this study we vary simultaneously factorisation and renormalisation scales, within 
a range of 0.5 to 2 times the reference central scale /i, where \i is the average of the two sparticle masses in the final 
state. In the future, we shall study the effect of varying renormalisation and factorisation scales independently. 
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The corrections often increase the size of the cross section with respect to the leading-order prediction [|8] |9] 
[Toll if the renormalisation and factorisation scales are chosen close to the average mass of the produced SUSY 
particles. As a result, the SUS Y-QCD corrections have a substantial impact on the determination of mass exclusion 
limits and would lead to a significant reduction of uncertainties on SUSY mass or parameter values in the case 
of discovery, see e.g. ifTTl . The processes listed in Eqs. ([TJ and (01 have been known for quite some time at 
NLO in SUS Y-QCD H2] [TJJ fH [JJl . Note that SUS Y-QCD corrections can be split into two parts: first, the 
QCD corrections induced by gluon or quark radiation and by gluon loops, which follow essentially the same 
pattern as, for example, the top pair production, and second, the virtual diagrams which involve squark and gluino 
loops, which are independent of the real emission corrections. For relatively heavy SUSY spectra the latter are 
numerically sub-leading terms which are, however, challenging to compute. This is mainly due to a large number 
of Feynman diagrams with many different mass scales contributing to the overall cross section. For example 
for stop pair production, where neither light-flavour squarks nor gluinos appear in the tree-level diagrams, these 
contributions can be easily decoupled lTT5l . The only part which requires some attention is the appropriate treatment 
of the counter term and the running of the strong coupling constant. This decoupling limit is implemented in 
Prospino2 ifTol . For light-flavour squark and gluino production this decoupling would only be consistent if 
applied to the leading order as well as NLO contributions. This is usually not required, unless we choose specific 
simplified models. 

Given the expected squark flavour structure in the MSSM, most numerical implementations, including PROSPIN02, 
make assumptions about the squark mass spectrum. The left-handed and right-handed squarks of the five light 
flavours are assumed to be mass degenerate. Only the two stop masses are kept separate in the NLO computations 
of light-flavour production rates ifPH [T_3] [14]. In the PROSPIN02 [ 16] implementation, this degeneracy is not as- 
sumed for the leading-order results. However, the approximate NLO rates are computed from the exact leading or- 
der cross sections times the mass degenerate K factors. For the pair production of third-generation squarks the four 
light squark flavours are mass degenerate, while the third generation masses are kept separate lfj"5l . This approxi- 
mation can for example be tested using MadGolem 07], an automized NLO tool linked to MadGraph4 |[T8l . 
It is also important to point out here that in PROSPIN02 the pair production of third-generation squarks is available 
as individual processes. However, sbottom pairs are included in the implicit sum of light-flavour squarks because 
there is no perfect separation of bottom and light-flavour decay jets. 

When summing the squark and gluino production rates including next-to-leading order corrections it is crucial to 
avoid double counting of processes. For example, squark pair production includes 0(af) processes of the kind 
qg —> qq*q. The same final state can be produced in qg production when the on-shell gluino decays into an anti- 
squark and a quark. The PROSPINO scheme for the separation and subtraction of on-shell divergences from the 
qq* process uniquely ensures a consistent and point-by-point separation over the entire phase space. For a finite 
particle mass it has recently been adopted by MC @NLO [ 19 1 for top quark processes. It is automized as part of 
MadGolem IfTTl . 

A significant part of the NLO QCD corrections can be attributed to the threshold region, where the partonic centre- 
of-mass energy is close to the kinematic production threshold. In this case the NLO corrections are typically large, 
with the most significant contributions coming from soft-gluon emission off the coloured particles in the initial 
and final state. The contributions due to soft gluon emission can be consistently taken into account to all orders 
by means of threshold resummation. In this paper, we discuss results where resummation has been performed at 
next-to-leading logarithmic (NLL) accuracy ll20l |2T1 l22l |231 l24l . 

The step from NLO to NLO+NLL is achieved by calculating the NLL-resummed partonic cross section cr( NLL ) 
and then matching it to the NLO prediction, in order to retain the available information on other than soft-gluon 
contributions. The matching procedure takes the following form 



where the last term in the square brackets denotes the NLL resummed expression expanded to NLO. The symbol 
{m 2 } stands for all masses entering the calculations and /i is the common factorisation and renormalisation scale. 
The resummation is performed in the Mellin moment N space, with all Mellin-transformed quantities indicated by 




+ 




(NLO) 



(3) 
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tilde. In particular, the Mellin moments of the partonic cross sections are defined as 

cry-Ma (N, {m 2 }, p 2 ) = [ dp p N - 1 aij^kl (P, {m 2 }, p 2 ) . (4) 
Jo 

The variable p = {mu + mi) 2 / s measures the closeness to the partonic production threshold and is related to 
the corresponding hadronic variable p — pXiXj in Eq. (O, where xi (xj) are the usual longitudinal momentum 
fraction of the incoming proton earned out by the parton The necessary inverse Mellin transform in Eq. (|3) 
is performed along the contour CT according to the so-called "minimal prescription" [1251 . The NLL resummed 
cross section in Eq. (f3]l reads 

Af LL > (N + l,Q 2 , p 2 ) Af LL > (N + 1, Q 2 , p 2 ) Ag^ L ] (N + l,Q 2 ,p 2 ) , (5) 

where the hard scale Q 2 is taken as Q 2 = (m^ + to/) 2 and d^j_^ kl / are the colour-decomposed leading-order 
cross sections in Mellin-moment space, with / labelling the possible colour structures. The functions A^ NLL ' ) and 
A^nll) sum the effects of the (soft-)collinear radiation from the incoming partons. They are process-independent 
and do not depend on the colour structures. These functions contain both the leading logarithmic as well as part of 
the sub-leading logarithmic behaviour. The expressions forAf LL) andAf LL) can be found in the literature [21 1. 
In order to perform resummation at NLL accuracy, one also has to take into account soft-gluon contributions 
involving emissions from the final state, depending on the colour structures in which the final state SUSY particle 
pairs can be produced. They are summarised by the factor 

A ( s ,NLL) (iV;Q2 ^ 2) = 

The one-loop coefficients Dj follow from the threshold limit of the one-loop soft anomalous-dimension matrix 
and can be found in |2"T1|221 . 

The analytic results for the NLL part of the cross sections have been implemented into a numerical code. The 
results of this code, added to the NLO results obtained from PROSPIN02, correspond to the matched NLO+NLL 
cross sections. Their central values, the scale uncertainty and the 68% C.L. pdf and a s uncertainties obtained using 
CTEQ6 . 6 J26) and MSTW2008 (27l PDFs have been tabulated for the squark and gluinoproduction processes 
of interest in the range of input masses appropriate for the current experimental analysisl£j The set of tabulated 
values has enough granularity in order to minimise the interpolation uncertainty when an intermediate value is 
used as input in the analysis. Together with a corresponding interpolation code, the tabulated values constitute the 
NLL-FAST numerical package 112411291 . 

Note that for specific processes, results beyond NLL accuracy are already available. The production of squark- 
antisquark pairs has been calculated at next-to-next-to-leading-logarithmic (NNLL) level 1351 . Also, a general 
formalism has been developed in the framework of effective field theories which allows for the resummation of 
soft and Coulomb gluons in the production of coloured sparticles [31 32] and subsequently applied to squark and 
gluino production 0211331 . In addition, the dominant next-to-next-to-leading order (NNLO) corrections, including 
those coming from the resummed cross section at next-to-next-to-leading-logarithmic (NNLL) level, have been 
calculated for squark-antisquark pair-production [34. 35 ]. The production of gluino bound states as well as bound- 
state effects in gluino-pair and squark-gluino production has also been studied ll36l [371 [3H1 [39l . Furthermore, 
electroweak corrections to the 0(a 2 ) tree-level processes [40. [41] |42] |43] 04] |45] [46] and the electroweak Bom 
production channels of 0(aa s ) and 0(a 2 ) 11471 l48l are in general significant for the pair production of SU(2)- 
doublet squarks q~L and at large invariant masses, but they are moderate for inclusive cross sections and will not be 
included in the results presented here. 

" Note that we use NLO pdfs with the NLO+NLL matched cross section calculation. The reduction of the factorisation scale 
dependence observed in the NLO+NLL predictions is a result of a better compensation between the scale dependence of the 
NLO evolution of the pdf and the short distance cross section, and does not depend on whether pdfs are fitted using NLO 
or NLO+NLL theory, see for example Ref. [28 1 . In general, one can apply NLL threshold resummation with NLO pdfs to 
processes like heavy SUSY particle production for which the summation of logarithms is more important than for the input 
data to the NLO fits. However, it would be interesting to systematically study the difference between NLO and NLO+NLL 
input to global pdf determinations for SUSY particle production at the LHC. 
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Resummation has also been studied for the production of electroweak SUSY partcles, see Refs. [49 50l l5"Tll52~ll53] 
|54j [55] [56) . However, since the electroweak processes do not significantly contribute to the inclusive SUSY cross 
section in general, we use the NLO-SUSY QCD calculation [57 58 59 1 for their evaluation. 

3 Treatment of cross sections and their associated uncertainties 

Here, the determination of the cross section central values and associated systematic uncertainties is discussed in 
detailed, in view of the recent progress discussed in the previous section. 

In case of pair production processes of strongly interacting particles, for which NLL calculations exist, the cross 
sections are taken at the next-to-leading order in the strong coupling constant, including the resummation of soft 
gluon emission at the NLL level of accuracy, performed using the NLL-FAST code. 

These cross sections are currently available for masses spanning 200 GeV to 2 TeV for squarks and gluinos and 
100 GeV to 1 TeV for direct stop or sbottom pair production. In the case of associated squark-gluino production 
and gluino-pair production processes, the calculations are extended up to squark masses of 3.5 TeV and 4.5 TeV, 
respectively. Following the convention used in PROSPIN02, in the case of squarks, which can be more or less 
degenerate depending on a specific SUSY scenario, the input mass used is the result of averaging only the first 
and second generation squark masses. Further details on different scenarios considered to interpret the variety of 
experimental searches developed by the ATLAS and CMS collaborations are described in Section IXTl 

In the case of other types of production processes, such as electroweak production and strong-weak production, or 
in few cases where the sparticle masses fall outside the range previously described, the NLO approach is considered 
using Prospino2. 

A special treatment is performed in scenarios in which either the squarks or the gluinos are set at a very large scale 
such that their production is not possible at the LHC. Defining this large scale is arbitrary and in some cases it may 
have a non-negligible impact in the production of the SUSY particles residing at the TeV scale (e.g. squarks at 
high scales can still contribute to the gluino pair production process via a i-channel exchange). Thus, instead of 
setting an explicit arbitrary value for the mass of the decoupled particle, the NLO+NLL calculation implemented 
in NLL-FAST assumes that this production does not interfere in any possible way with the production processes of 
the rest of the particles. Thus, the particle is completely decoupled from the rest of the phenomenology at the TeV 
scale. 

The uncertainties due to the choice of the renormalisation and factorisation scales as well as the PDF are obtained 
using the NLL-FAST code or computed with PROSPIN02. In order to combine all these predictions and obtain 
an overall uncertainty, the PDF4LHC recommendations are followed as close as possible, subject to the current 
availability of the different calculations. Thus, an envelope of cross section predictions is defined using the 68% 
CL. ranges of the CTEQ6 . 6 J26) (including the a s uncertainty) and MSTW2008 [27] PDF sets, together with 
the variations of the scales. The nominal cross section is obtained using the midpoint of the envelope and the 
uncertainty assigned is half the full width of the envelope. Mathematically, if PDFcteq,up (PDFcTEQ,down) 
and SCActeq,up (SCAcTEQ,down) are the upwards (downwards) one sigma variations of the CTEQ6 . 6 PDF 
set, respectively, PDF M stw,u P (PDF M STW,down) and SCA M stw,u P (SCA M STW,down) are the corresponding 
variations for the MSTW2 008 PDF set and, finally, as, U p (as, down) is the corresponding up (down) one sigma 



uncertainty of the a s coupling constant, the following quantities can be calculated: 

CTEQ up = ^/PDF2 TEQ up + SCA^ TEQiUp + a|_ up , (7a) 

CTEQ down = ^/PDFcTEQ,down + SCAcTEQ,down + a| down , (7b) 

MSTW up - ^/PDFgjsTw^p + SCA^Tw^p , (7c) 

MSTWdown = ^/PDF^gTW.down + SCA^STW.down • (7d) 

The corresponding upper and lower values of the envelope created by this set of numbers and the nominal predic- 
tions (CTEQ nom and MSTW nom ) is obtained by: 

U = max(CTEQ nom + CTEQ up , MSTW nom + MSTW up ) , (8a) 

L = min(CTEQ nom - CTEQ down , MSTW nom - MSTW down ) , (8b) 
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and the final corresponding cross section (ct) and its symmetric uncertainty (Act) is taken to be: 



Act 



(U + L)/2, 
(U-L)/2. 



(9a) 
(9b) 



Full compliance with the PDF4LHC recommendations, with the inclusion of other PDF sets such as NNPDF |60], 
will be implemented in the near future. We notice that, as discussed in section|4] the additional contribution to the 
systematics coming from a s uncertainties is negligible. 

3.1 Special cases 

Some models require special treatment in order to make sure that the NLO cross sections are correctly computed. 
This is important because of the way higher order calculations handle different squark flavours, which could lead 
to double-counting of diagrams. Given the difficulty to provide a comprehensive summary of all situations that are 
being considered in the interpretation of the LHC data, we only discuss here few relevant cases, to exemplify the 
approach followed. 

Simplified Models 

A variety of simplified models [61 1 are considered by the experiments. In some cases, the gluino, sbottom, and 
stop squarks are decoupled from the rest of the supersymmetric spectrum^ In this specific simplified model, only 
squark-antisquark production is allowed and this process is flavour-blind, if the masses are considered degenerate. 
Since the NLO+NLL calculations consider the sbottom as degenerate in mass with the squarks of the first and 
second generations, the overall cross section has to be rescaled down by a factor of 4/5 = 0.8. 

In other cases where the gluino is not decoupled, squark-gluino and squark-squark productions are feasible and they 
are not corrected further. The only effect could come from a fe-quark in the initial state, which is suppressed [23 1 . 
Other types of simplified models decouple not only the squarks from the third generation particles, but also all 
the right-handed squarks. These scenarios primarily focus on final state decays via charginos or neutralinos. The 
squark mass is calculated by averaging the non-decoupled squark masses and the final cross section is rescaled 
down by a factor of (4/5) • (1/2) = 0.4. 

Treatment of 3rd generation squarks 

Direct stop and sbottom production must be treated differently from the rest of squark families because, for in- 
stance, the i-channel gluino-exchange diagrams are suppressed. In beyond the leading-order computations pro- 
cesses involving sbottoms can be degenerate with the rest of the other squark flavours, this could lead to a potential 
double counting of sub-processes. In order to avoid this, in scenarios in which the production of different squark 
flavours are present, the squark pair production cross section is rescaled down to subtract the sbottom contribution 
and the corresponding process is computed separately. 

At leading order the production cross sections for the third -generation squarks depend only on their masses, and 
the results for sbottom and stop of the same mass are therefore equal. At NLO in SUSY-QCD, additional SUSY 
parameters like squark and gluino masses or the stop/sbottom mixing angle enter. Their numerical impact, however, 
is very small JT5JH3J. A further difference between stop and sbottom pair production arises from the bb — > bb* 
channel, where the initial-state bottom quarks do allow a t-channel gluino-exchange graph that gives rise to extra 
contributions. However, as has been demonstrated in Ref. Il23l their numerical impact on the hadronic cross 
sections is negligible. Thus, for all practical purposes, the LO and higher-order cross-section predictions obtained 
for stop-pair production apply also to sbottom-pair production if the input parameters, i.e. masses and mixing 
angles, are modified accordingly. 



The production cross sections and associated uncertainties resulting from the procedure described in the previous 
section are discussed here for different processes of interest. We discuss three distinct cases: gluino pair production 

2 ' While this scenario is possible if we only consider the TeV scale, it bears some challenges at higher energy scales. Any kind 
of renormalisation scale evolution will generate squark masses at the scale of the gluino mass, but not vice versa. Thus, any 
discovery of light squarks associated with heavy gluinos would point to a non-standard underlying model [62 1 . In spite of 
all theory prejudice it is clearly adequate that these regions be experimentally explored. 



4 Gluino and squark production at the LHC 
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with squarks decoupled, squark-antisquark pair production with gluino decoupled and stop/sbottom pair produc- 
tion. The results shown here are mainly illustrative: tables with cross sections and systematics obtained in other 
scenarios which take into account the full complexity of the spectrum, as a function of the parameters of various 
SUSY models, are collected at the SUSY cross section working group twiki page 0. 

4.1 Gluino pair production 

The gluino pair production cross section in a model where the squarks are decoupled is shown in Figure Q] for a 
gluino mass range spanning the current sensitivity of the ATLAS and CMS experiments. In the figure, the black 
(red) line corresponds to the NLO+NLL nominal cross section and renormalisation and factorisation scale uncer- 
tainties obtained using the CTEQ6 . 6 (MSTW2 8) PDF set. The solid yellow (dashed black) band corresponds to 
the total uncertainty of the cross section using CTEQ6 . 6 (MSTW2 8), as derived from Eq.|7a] Finally, the green 
solid lines delimit the envelope and the central value. They correspond to the central nominal value along with the 
total uncertainties. 




Figure 1: NLO+NLL gluino pair production cross section with squarks decoupled as a function of mass. The dif- 
ferent styled black (red) lines correspond to the cross section and scale uncertainties predicted using the CTEQ6 . 6 
(MSTW2008) PDF set. The yellow (dashed black) band corresponds to the total CTEQ6 . 6 (MSTW2008) uncer- 
tainty, as described in the text. The green lines show the final cross section and its total uncertainty. 

As it can be seen, the cross sections predicted using the MSTW2 008 PDF set are lower than the ones predicted by 
CTEQ6 . 6. For this mass range the uncertainty using MSTW2008 is approximately half that of the CTEQ6 . 6 set. 
Thus, the impact of the MSTW2 8 prediction in the current prescription is almost negligible. The small difference 
between the green and black dotted lines also shows that the impact of the a s uncertainty in the CTEQ6 . 6 case 
is negligible. This is consistent with the intuition that at the large values of x that are relevant for SUSY particle 
production, the intrinsic PDF uncertainty dominates over the small evolution differences caused at these high scales 
by the a s uncertainty. This justifies as a good approximation not including the a s uncertainty in the MSTW2 00 8 
case. 

4.2 Squark-antisquark production 

In order to show the evolution of the squark-antisquark production cross section as a function of the squark mass, 
a scenario has been chosen in which the gluino is decoupled. The results are shown in Figure |2] using the same 
graphical convention described in the gluino pair production case. The central values predicted using the CTEQ6 . 6 
and MSTW2 8 PDF sets are relatively close to each other. Again, the MSTW2 8 uncertainty is the smallest of 
the two, although relatively larger than in case of gluino pair production. 
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Figure 2: NLO+NLL squark-antisquark production cross section with gluinos decoupled as a function of mass. 
The different styled black (red) lines correspond to the cross section and scale uncertainties predicted using 
the CTEQ6.6 (MSTW2 8) PDF set. The yellow (dashed black) band corresponds to the total CTEQ6.6 
(MSTW2 008) uncertainty, as described in the text. The green lines show the final cross section and its total 
uncertainty. 

4.3 Direct stop and sbottom productions 

The production cross section as a function of the stop mass for a model in which only the lightest stop is reachable 
by the LHC is shown in Figure [3] It should be noted that these cross sections are the same to those of a model in 
which only the lightest sbottom is accessible, assuming the rest of the coloured SUSY spectrum decoupled. 

We discuss two different cases depending on the stop mass. For low stop masses, the cross sections predictions 
using the MSTW2 008 PDF set are larger than the ones using CTEQ6 . 6, but the relative uncertainties are similar. 
Thus, the total cross section and the overall uncertainty are larger than those predicted by CTEQ6 . 6 alone. For 
larger stop mass ranges, the predictions using any of the two PDF sets under consideration are quite similar for the 
central value and its associated uncertainties. The final prediction is mostly determined by the CTEQ6 . 6 PDF set. 

5 SUSY electroweak production at the LHC 

SUSY electroweak particles are usually produced via cascade decays but they can also be directly produced. In 
this section we study the production cross sections along with the theoretical uncertainties using a CMSSM ll63l 
scenario with tan/3 = 10, Aq = GeV and /i > 0. Figure [4] shows the NLO cross section for various SUSY 
electroweak productions in the mo — mi/a mass plane. The pure electroweak states have much lower cross 
sections in comparison with SUSY colour production. All possible combinations of charginos/neutralinos (x) and 
sleptons (J) are considered for the rate given by \X an d H sub-processes, respectively. The theoretical uncertainty 
due to the scale variation is found to be less than 10% for pure electroweak states as shown in Figure [5] The 
squark and electroweak gaugino associated sub-processes can have errors as large as 35%. The systematics due 
to the CTEQ6 . 6 and MSTW2008 PDFs for XX are found to be less than 5%. The uncertainties on the rest of the 
sub-processes, as shown in Figure|6]and Figure |7] are negligible. 

6 Summary and future prospects 

We presented here a proposal^ based on state-of-the-art higher-order calculations, for the definition of benchmark 
production cross sections and associated theoretical uncertainties of various SUSY processes of interest at the 

3 ' This documentation should be cited along with the original papers, i.e. Refs. [14 20 .21, 22 24 1 for inclusive squark/gluino 
production, Refs. [ 15 23 24] for stop or sbottom direct production, and Ref. 1 57 1 for the production of electroweak gauginos. 
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Figure 3: NLO+NLL stop-antistop production cross section as a function of mass. The different styled black (red) 
lines correspond to the cross section and scale uncertainties predicted using the CTEQ6 . 6 (MSTW2 08) PDF set. 
The yellow (dashed black) band corresponds to the total CTEQ6 . 6 (MSTW2008) uncertainty, as described in the 
text. The green lines show the final cross section and its total uncertainty. The figure on the top (bottom) shows a 
low (high) stop mass range. 
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Figure 4: NLO cross sections for SUSY electroweak productions in the mo — tuxii plane using CTEQ6 . 6 PDF 
set. The CMSSM framework is used, with tan (3 = 10, A = GeV and fx > 0. 
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Figure 5: Uncertainties in percentage due to the scale choice for SUSY electroweak productions in the mo — mi/2 
plane using CTEQ6 . 6 PDF set. The CMSSM framework is used, with tan /3 = 10, A = GeV and [i > 0. 
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Figure 6: Uncertainties in percentage due to CTEQ6 . 6 PDF error prescription for SUSY electroweak productions 
in the too -mi/j plane. The CMSSM framework is used with tan (3 = 10, Aq — GeV and /j, > 0. 
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Figure 7: Uncertainties in % due to MSTW2 008 PDF error prescription for SUSY electroweak productions in the 
mo — mi/2 plane. The CMSSM framework is used, with tan f3 = 10, Aq = GeV and [i > 0. 
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LHC. The common set of prescriptions described in this document aims to establish a reference framework for the 
theoretical input used by the ATLAS and CMS collaborations in the interpretation of their measurements. 

We presented a collection of explicit results for collisions at a centre-of-mass energy of 7 TeV. The theoretical 
systematic uncertainties are larger for higher masses, and they are typically dominated by the PDF uncertainties. 
These have a significant impact when assessing the experimental constraints or the sensitivity to a given SUSY 
model. Detailed numbers and tables for a broad class of SUSY models and parameters are collected at the SUSY 
cross section working group twiki page [7 ], where the results corresponding to higher collision energies will also 
appear soon. 

The LHC has recently raised the centre-of-mass energy and the experiments are already analysing data from these 
collisions. Hence, the current set of prescriptions will evolve and is expected to adapt to 8 TeV collisions. In 
particular, the mass ranges for which the NLO+NLL calculations are available will be enlarged including usage of 
the complete PDF4LHC prescription. 
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